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We discuss the construction and properties of an approximate solution of the Ginsparg- Wilson equation, the 
so-called chirally improved lattice Dirac operator. In particular we study the behavior of its eigenmodes in smooth 
instanton backgrounds as well as for thermalized gauge configurations on both sides of the QCD phase transition. 
We compare with results from other Dirac operators including the overlap operator. The results support the 
picture of chiral symmetry breaking being closely related to instantons. 



1. Motivation and introduction 

During the last three years it has been real- 
ized that the Ginsparg- Wilson relation (IJ is the 
key to chirally symmetric fermions on the lattice. 
However, the overlap operator which exactly 
solves the Ginsparg- Wilson equation is numeri- 
cally very expensive. For many applications an 
alternative to an exact solution is an approximate 
Ginsparg- Wilson Dirac operator, such as a finite 
parametrization of the perfect action [^| or the re- 
cently proposed chirally improved Dirac operator 
0. The latter results from a systematic expan- 
sion of a solution of the Ginsparg- Wilson equa- 
tion. 

In this contribution we report on our results 
for the eigenvectors and eigenvalues of the chirally 
improved Dirac operator, using smooth instanton 
backgrounds as well as thermalized (quenched) 
gauge configurations on both sides of the QCD 
phase transition. These studies serve two pur- 
poses: Firstly, they allow to obtain a sound un- 
derstanding of the properties of the chirally im- 
proved Dirac operator. In particular we also com- 
pare our results to calculations done with the 
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overlap operator in order to understand the ef- 
fects of small violations of chirality. Secondly, the 
analysis of the eigenmodes provides insights into 
the nature of the QCD vacuum, since the eigen- 
vectors couple to localized objects in the under- 
lying gauge field such as instantons or instanton 
anti-instanton molecules. 

2. Remarks on the construction of the lat- 
tice Dirac operators 

Chirally improved operator: In Q we sug- 
gested an expansion of the lattice Dirac operator 
D in terms of a systematic series of simple hop- 
ping terms, connecting sites along paths of arbi- 
trary length and allowing for all elements of the 
Clifford algebra. For the ansatz one respects all 
standard translational and rotational symmetries 
of the lattice Dirac operator as well as invariance 
under C and P and 75-hermiticity. We insert that 
formal series into 

E= -D - D f + D ] D . (I) 

Finding a solution to the Ginsparg- Wilson equa- 
tion jjj amounts to solving for E = 0. The prod- 
uct D' D is again constructed algebraically, lead- 
ing to terms involving Clifford algebra elements 
and paths corresponding to ordered products of 
gauge link variables. 
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E is an infinite series of terms that are linear 
and quadratic in the coupling coefficients of the 
Dirac operator multiplying the path terms. In- 
dependent terms have to vanish and thus we ob- 
tain a system of coupled equations quadratic in 
the (unknown) coupling constants. We truncate 
the ansatz for D and the series for E at paths 
of length 4. Whereas the complete series would 
lead to an exact solution of the Ginsparg- Wilson 
equation, the solution now will only be approxi- 
mate and it remains to be studied, how good this 
approximation is. 

In solving the set of equations there is much 
freedom. In order to partially remedy the approx- 
imation we introduce two more conditions (and 
parameters) such that the boundary condition for 
a massless Dirac operator 

D(p) = ii/+0{ P 2 ) , (2) 

is satisfied at a given value of the gauge coupling. 
The solution we study here has 19 terms and val- 
ues of the coefficients are listed in the appendix of 
||. It is ultralocal with non- vanishing couplings 
only for \x — y\ < y/5 decreasing exponentially in 
size. 

Exact solutions to the Ginsparg- Wilson equa- 
tion (E = 0) have eigenvalues on a unit circle with 
center 1 in the complex plane. It turned out that 
our solution, defining a chirally improved Dirac 
operator, has eigenvalues, that are very close to 
the unit circle in particular for the small eigen- 
value region, which is most relevant for the in- 
frared behavior. 

Overlap operator: Below we compare the 
properties of the sector with small eigenvalues 
(complex or real) of the introduced chirally im- 
proved operator with those of an overlap operator 
|2j which has the form 

D ov =l-Z with Z = 75 sign(tf) , (3) 

where H is related to the hermitian Dirac opera- 
tor. This is constructed from an arbitrary Dirac 
operator D (e.g. the Wilson operator), 

H = 75 (a - D Q ) , (4) 

and s is a parameter which may be adjusted in 
order to minimize the probability for zero modes 



of H. In case Dq is already an overlap opera- 
tor one reproduces D ov = Dq for s = 1 because 
sign (sign(iJ)) = sign(iJ). 

The sign-function may be defined through the 
spectral representation 

sign H = ^2 si 8' n ^ I*) (* I > ( 5 ) 

i 

(|z) denote the eigenvectors) but in practical com- 
putations this definition cannot be used for real- 
istic QCD Dirac operators, which for L 4 lattices 
have dimension O(10 5 — 10 6 ) (~ n co ior ' nDirac ■ 
L 4 ). Note that in the subsequent applications the 
operator has to be applied many times, since it 
may be entering a diagonalization or a conjugate 
gradient inversion tool. One therefore relies on 
the relation 



sign(ff) = 




and approximates the inverse square by some 
method. In our computations |6) we follow the 
methods discussed in j^g • One approximates the 
inverse square root by a Chebychev polynomial, 
which has exponential convergence in [e, 1], where 
e (and thus the order of the polynomial) depends 
on the ratio of smallest to largest eigenvalue of 
H 2 . Clenshaw's recurrence formula provides fur- 
ther numerical stability. 

Depending on the input Dirac operator Dq the 
rate of convergence may be quite unfavorable. Al- 
though one may adjust s it turns out that partic- 
ular simple operators like the Wilson Dirac opera- 
tor eventually will give rise to several small eigen- 
values of H 2 . Technically, one then proceeds as 
follows: One determines the subspace of e.g. the 
lowest 20 eigenmodes of H 2 and computes the in- 
verse square separately for the subspace (using 
the spectral representation) and the reduced op- 
erator (with polynomial approximation) . For this 
approach it is important to determine the sub- 
space with high accuracy (we request 12 digits). 

In our study we need to diagonalize not only 
the hermitian matrix H but, for the later analysis, 
also D ov and the chirally improved operator. All 
diagonalizations have been done with the Arnoldi 
method ||. 

CPU-time comparison: The central element 
of the diagonalization (and conjugate gradient 
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methods) is the multiplication of the Dirac op- 
erator with a vector. For the consideration of 
the computational effort we disregard the startup 
time for initialization of the operator. Comparing 
the CPU-time in units of the time used for the 
Wilson Dirac operator we find that the chirally 
improved operator needs 24 units and the over- 
lap operator (with Wilson input) typically needs 
300±100 units (always with 12 digits accuracy for 
the matrix norm | \H 2 — 1| |). These are the results 
for lattice size 12 4 , but they are not significantly 
different for 8 4 or 16 4 . 

As mentioned, whenever (1 — D$) is already an 
overlap operator, one finds D ov — D - We are in- 
terested in the case when (1 — D ) is close to, but 
not quite an overlap operator. Such an operator 
might accelerate the convergence in the construc- 
tion of the overlap operator, since the conver- 
gence of the polynomial approximation depends 
crucially on the smallest eigenvalue of H 2 . It is 
therefore of high practical importance, whether H 
will have eigenvalues close to zero. A Dirac op- 
erator D used as input to the construction of H 
may produce such small eigenvalues. Physically 
they are related to the appearance or disappear- 
ance of an instanton. Different Dirac operators 
will differ in their sensitivity to "identify" such 
structures in the gauge configuration. 

When using the chirally improved operator - 
in comparison to using the Wilson Dirac opera- 
tor - as starting point we need typically a factor 2 
less matrix- vector multiplications (75 ± 25 terms 
of the Chebychev series for 12 4 or 16 4 lattices). 
One reason that this improvement factor is not 
higher may lie in the removal of the small eigen- 
mode subspace, which is relatively more efficient 
for the "bad" Wilson operator used as input to 
the overlap construction. Another reason may be 
that the series approximation forces all eigenval- 
ues to the unit circle; since the chirally improved 
operator has its large eigenvalues not as close to 
the circle as the small ones, one has to pay the 
price with longer polynomials. 

3. Small eigenvalues and instantons 

For an instanton in the continuum it is known 
that the Dirac operator has a zero mode. On the 
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Figure 1. Position x of the real eigenvalue (the 
would-be zero mode) of different lattice Dirac op- 
erators as a function of the radius R of the un- 
derlying instanton configuration. 



lattice only an exact solution of the Ginsparg- 
Wilson equation will display an eigenvalue which 
vanishes exactly. An approximate Ginsparg- 
Wilson Dirac operator will then have a small real 
eigenvalue. Typically the size of the real part 
will increase as one shrinks the radius of the un- 
derlying instanton. The rate at which the real 
eigenvalue moves away from the origin is a mea- 
sure for the quality of the approximate solution of 
the Ginsparg- Wilson equation. A good approx- 
imation will keep the eigenvalue near zero also 
for relatively small instantons. In order to im- 
plement this test of the would-be zero modes we 
use instantons with several values of the radius 
R discretized on the lattice following the method 
described in |ll]||. 

Fig. shows the position x of the real eigen- 
value (the would-be zero mode) of different Dirac 
operators as a function of the radius R of the 
underlying instanton. We show the results com- 
puted on 16 4 lattices for the chirally improved 
operator (circles), the Wilson operator (squares), 
and the overlap operator (diamonds). The lat- 
ter, since it is an exact solution of the Ginsparg- 
Wilson equation, has an exact zero mode. It is, 
however, interesting to see that the overlap oper- 
ator can trace the instanton only down to a radius 
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Figure 2. Probability distribution of the position 
x of the real eigenvalues of the chirally improved 
Dirac operator. The data were computed on 16 4 
lattices at values of j3\ = 8.10, 8.30 and 8.45 using 
200 configurations for each j3\. 



R = 1.5. The Wilson operator has an eigenvalue 
which very quickly departs from zero as R is de- 
creased. The chirally improved operator keeps 
the real eigenvalue fairly close to even for rel- 
atively small instantons and only for R < 2 the 
deviation becomes significant. 

For practical applications of the chirally im- 
proved Dirac operator, the question is whether 
also in a gauge configuration with quantum fluc- 
tuations the real eigenmodes lie close to the ori- 
gin. In order to test this property, we computed 
histograms for the probability distribution of the 
real eigenvalues. The gauge fields were gener- 
ated in a quenched simulation, using the tadpole 
improved Luscher-Weisz action pTlJ . The his- 
tograms were computed on 16 4 lattices at values 
of Pi = 8.10, 8.30 and 8.45 using 200 configura- 
tions for each [3\. The corresponding values of 
the couplings 02 and /?3 of the gauge action can 
be found in p2[ . There we obtained for the lat- 
tice spacing a the values 0.127 fm, 0.107 fm, and 
0.1 fm for [3i = 8.10, 8.30, and 8.45, respectively. 

From Fig. ^ it is obvious that the eigenvalues 
are concentrated near the origin. The tail of the 
distribution decays very quickly and extends al- 
most exclusively towards values x > 0. Both 



these features are highly welcome: The narrow 
width of the distribution lets one hope for only a 
very small additive renormalization of the quark 
mass. The fact that the distribution does not ex- 
tend to values x < implies furthermore that 
exceptional configurations (which spoil the com- 
putation of the propagator) are suppressed effec- 
tively. 

4. Localization properties of zero modes 

We now analyze the eigenvectors with real 
eigenvalues, i.e. the modes which correspond to 
the zero modes in the continuum. For a contin- 
uum instanton the zero modes are known to be 
localized at the same position as the instanton. 
This property also holds for thermalized configu- 
rations on the lattice [pUf . The eigenmodes pro- 
vide an efficient filter for the excitations of the 
QCD vacuum as seen by the Dirac operator. 

To study the eigenvectors we define the scalar 
density po(x). Let il>(x) c be an eigenvector of the 
Dirac operator with x denoting the lattice point 
while c is the color index. A gauge invariant den- 
sity p a {x) is defined as 



ip(x)* c Y a ij;(x) c 



(7) 



(Here Tq = 1, Tg =75; the sum over Dirac indices 
is implied.) In the last section we have seen that 
the continuum zero modes correspond to eigen- 
vectors of the lattice Dirac operator with eigen- 
value zero (overlap operator) or small real eigen- 
value (chirally improved operator, Wilson opera- 
tor). An interesting question is, to which extent 
the eigenvectors are different for different lattice 
Dirac operators. 

In Fig. H we compare iso-surfaces of the scalar 
density for the zero-mode of different Dirac oper- 
ators on the same gauge configuration (12 4 lattice 
at Pi = 8.45). All three operators see the same 
localized instanton. Only the details vary slightly 
as can be seen by the different fluctuations around 
the core. 

An interesting effect was observed for small in- 
stantons in Q where we compared the localiza- 
tion of the zero mode for the smooth instantons 
discussed above for different lattice Dirac opera- 
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Figure 3. Comparison of iso-surfaces (for a given 
time-slice) of the scalar density for the zero-mode 
of different Dirac operators. We show the den- 
sity for the Wilson operator (top), the chirally 
improved operator (middle), and the overlap op- 
erator (bottom). 
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Figure 4. Comparing the localization of the zero 
modes of the chirally improved operator and the 
overlap operator. 



tors. It was found that for instantons with radius 
smaller than 3 in lattice units, the overlap oper- 
ator has a zero mode which is less localized than 
in the case of a continuum instanton. Since the 
other two operators (chirally improved, Wilson) 
do not display such a blowing up of zero modes we 
attribute this effect to the rather big size of the 
overlap operator. As opposed to the other two 
operators the overlap operator is not ultra-local, 
i.e. \D ov (x,y)\ decays exponentially with \x — y\ 
but does not vanish exactly. Objects smaller than 
the typical size of the overlap operator thus can 
no longer be resolved. We estimate j(| that at 
typical lattice spacings of a ~ 0.1 fm this effect 
sets in at about 0.3 fm. 

For further study we introduce the so-called in- 
verse participation ratio, a convenient measure of 
the localization widely used in solid state physics. 
Due to the normalization of the eigenmodes we 
have J2 X Po( x ) — 1- The inverse participation 
ratio is then defined as 

7 = v]Tp (z) 2 , (8) 

X 

where V is the volume of the lattice. An alterna- 
tive observable for localization, based on the self 
correlation of the scalar density po(x) was ana- 
lyzed in @. 
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In Fig. [| we show a scatter plot of the inverse 
participation ratio for zero modes of the overlap 
operator I ov plotted as a function of the inverse 
participation ratio of the corresponding modes of 
the chirally improved operator I c h. imp. ■ The data 
were computed on 12 4 lattices using 200 config- 
urations at ft = 8.10 (circles), and ft = 8.45 
(filled triangles). The full curve with filled circles 
is the result for the smooth artificial instantons. 
The dashed line is a straight line with slope 1. 
If both operators would see the same localization 
for the zero modes, all data points would be lo- 
cated on this line. However, obviously this is not 
the case. Except for a single data point, all data 
are below this line. This means that the over- 
lap operator systematically produces zero modes 
with smaller localization when compared to the 
chirally improved operator; this trend becomes 
stronger for more localized states. 

5. Analysis of the near-zero modes 

So far we have concentrated on analyzing the 
eigenvectors with real eigenvalues, i.e. the would- 
be zero modes. However, also the eigenvectors 
with small but complex eigenvalues, the so-called 
near-zero modes, play an important role in QCD 
phenomenology. In particular the eigenvalue den- 
sity p(X) near at the origin is related to the chi- 
ral condensate through the Banks-Casher formula 

in, 



(ipip) = — 7r lim lim p(X)V 1 . 



(9) 



According to the instanton liquid picture of chi- 
ral symmetry breaking p6| , the near-zero modes 
come from weakly interacting instantons and 
anti- instantons. The corresponding eigenvectors 
are expected to be still quite localized, since they 
correspond to interacting "lumps" which still re- 
semble instantons or anti-instantons. These pre- 
dictions may be tested in an ab-initio calculation 
on the lattice. Again we use the inverse partici- 
pation ratio defined in Eq. (|j) to study the local- 
ization of these modes. An interesting question 
is how the localization of eigenvectors varies as a 
function of the eigenvalue. 

In Fig. |5| we show a histogram of the inverse 
participation ratio for eigenvectors as a function 
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Figure 5. Histograms of the inverse participation 
ratio for eigenvectors as a function of the imagi- 
nary part of the corresponding eigenvalue. 



of the imaginary part of the corresponding eigen- 
values. The data were computed on 16 4 lattices 
at values of ft = 8.10,8.30 and 8.45 using 200 
configurations for each ft. It is obvious that the 
most localized states have eigenvalues near the 
origin. Such large values are expected for weakly 
interacting instantons and anti-instantons. As 
Im A is increased, the localization decreases. The 
corresponding modes do no longer feel distin- 
guished topological lumps but become dominated 
by quantum fluctuations. 

An interesting observable was suggested by 
Horvath et al. in jL7|. Their local chirality vari- 
able was subsequently studied by several other 
groups in jl8|,| 12 1. The underlying idea is to test 
whether the near-zero modes are locally chiral, 
i.e. left-handed for lattice points near the instan- 
ton peak of the density and right-handed near the 
anti-instanton peak. The construction starts with 
defining left- and right-handed densities p±(x), 

P±(x) = ^(po(x)±p 6 (x)) , (10) 



i.e. the density (Q) projected on positive and 
negative chirality. A mode dominated by in- 
stantons and anti-instantons is expected to have 
p+ (x) > and p_ (x) = for lattice points x near 
the instanton peak and vice versa for x near an 
anti-instanton peak. Thus the ratio p + (x)/p_(x) 
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Figure 6. Local chirality of eigenvectors as a func- 
tion of Im A of the corresponding eigenvalues. We 
use cuts of 1% (top) 6.25% (middle) and 12.5 % 
(bottom) for the number of supporting points. 



should be large for x close to an instanton peak 
and small for all x close to an anti-instanton peak. 
In a final step Horvath et al. map the positive real 
values of this ratio onto the interval between —1 
and +1 using the transformation 



X(x) = -arctan [J^T 

7T \y P-{X) 



1 . 



(11) 



In order to suppress the quantum fluctuations one 
evaluates X{x) only for x near instanton or anti- 
instanton peaks, i.e. for x near the maxima of 
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Figure 7. The inverse participation ratio as a 
function of the imaginary part of the correspond- 
ing eigenvalue. These data were computed at 
Pi = 8.45 where (with time extent 6) the the- 
ory is in the deconfined phase. In the subsample 
used here the expectation value of the Polyakov 
loop (L) is close to real. 



the scalar density po(x). We computed our re- 
sults using cutoffs of 1%, 6.25% and 12.5% for 
the number of lattice points x (ordered according 
to po{x)) and averaged over X{x). On the gauge 
field ensemble the local chirality variable is ex- 
pected to show a double peak structure only for 
instanton dominance. 

In Fig. ^| we show the local chirality of the near- 
zero modes as a function of the imaginary part of 
their eigenvalues. The data were computed from 
200 configurations on 16 4 lattices at (3\ = 8.45. 
When using a cut of 1% on the number of sup- 
porting lattice points (top plot) we find the best 
signal while for cuts of 6.25% (middle plot) and 
12.5% (bottom plot) the signal becomes weaker 
due to quantum fluctuations. Furthermore it is 
obvious for all three plots that the strongest sig- 
nal for local chirality is found for the eigenvectors 
with eigenvalues closest to the origin. For larger 
values of |ImA| the signal of local chirality for 
the corresponding eigenvectors weakens and the 
mode becomes dominated by quantum fluctua- 
tions. 
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6. The near-zero modes at finite tempera- 
ture 

Let us now discuss what happens to instan- 
tons at the QCD phase transition. At the crit- 
ical temperature T c the theory enters the decon- 
fincd phase and chiral symmetry is restored. The 
spectrum develops a gap (see e.g. Q for a plot 
of typical spectra) and the density of eigenvalues 
is zero near the origin. According to the Banks- 
Casher formula (^|) this gap in the spectral density 
amounts to a vanishing chiral condensate. 

The topological charge does not vanish 
abruptly at T c and topological objects are ob- 
served also for T > T c . Thus the liquid of weakly 
interacting instantons responsible for chiral sym- 
metry breaking at T < T c must undergo some 
change in its structure as the temperature is in- 
creased above T c . It may be expected that above 
T c instead of interacting weakly, instantons and 
anti- instantons form tightly bound molecules JTgJ . 
This prediction from instanton models can be 
tested on the lattice (compare [Jill for a similar 
study using the staggered Dirac operator). 

In Fig. |t] we show the inverse participation ra- 
tio as a function of the imaginary part of the 
corresponding eigenvalue, as in Fig. |, but now 
for finite temperature. In particular we show 
the results for configurations where the Polyakov 
loop (L) is in its real branch, which is the do- 
main of the Polyakov loop value for the full, un- 
qucnchcd theory. For the equivalent plot with 
complex Polyakov loop see ||. Due to the spec- 
tral gap in the deconfined phase, there are no 
near-zero modes and the graph has support for 
|Im A| > 0.15. One finds that now the most local- 
ized states are near the edge of the spectrum. Lo- 
calization decreases quickly as |ImA| is increased 
and the mode becomes a bulk mode dominated 
by quantum fluctuations. The inverse participa- 
tion ratio, however, provides only information on 
the localization of the modes and does not test for 
other properties such as local chirality. Thus it 
makes sense to repeat last section's study of Hor- 
vath et al.'s local chirality variable for T > T c . 

In Fig. H we compare the local chirality of 
the near-zero modes for an ensemble at T < T c 
(top plot, 6 x 20 3 lattices, 400 configurations at 
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Figure 8. Comparison of the local chirality of the 
(complex) near-zero modes for T <T C (top) and 
T>T C (bottom). 



/?! = 8.10) with an ensemble at T > T c (bot- 
tom plot, 6 x 20 3 lattices, 147 configurations with 
real branch Polyakov loop at (3\ — 8.45). The 
plots were computed using only the eigenvectors 
with |ImA| less than some cut corresponding to 
roughly the 4 eigenvalues closest to the origin, re- 
spectively the edge of the spectrum. Zero modes 
were left out. For T < T c we find a very pro- 
nounced double peak structure, i.e. a strong sig- 
nal for local chirality. At T > T c the double peak 
structure is gone and we do not observe local chi- 
rality in the eigenmodes. It is obvious that the 
weakly interacting instantons and anti-instantons 
undergo a drastic change as T increases above T c . 
The resulting states are still localized (see Fig. ^) 
but the local chirality present for weakly interact- 
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ing topological excitations has vanished, probably 
due to the formation of tightly bound molecules. 
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